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ABSTRACT 


Vinti has recently modified his spheroidal potential so 
as to permit exact inclusion of the effects of the third zonal 
harmonic of the planet's gravitational field. This cor- 
responds to a potential fitted exactly through the third zonal 
harmonic and about two-thirds of the fourth. The present 
paper treats the method for obtaining the position and ve- 
locity coordinates of a satellite moving in the field corre- 
sponding to this accurate modified potential. 
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COMPUTATIONAL PROCEDURE FOR VINTI’S 
ACCURATE REFERENCE ORBIT WITH INCLUSION 
OF THE THIRD ZONAL HARMONIC 


by 

N, L. Bonavito 
Goddard Space Flight Center 


INTRODUCTION 

Vinti (Reference 1) has found a gravitational potential for an axially symmetric planet in oblate 
spheroidal coordinates. This solution accounts for all of the second zonal harmonic and more than 
half of the fourth zonal harmonic. This potential, which simultaneously satisfies Laplace’s equa- 
tion and separates the Hamilton Jacobi equation, succeeds in reducing the problem of satellite 
motion to quadratures. 

More recently however, Vinti (Reference 2) has generalized his potential by means of a 
metric preserving transformation of the associated Cartesian system. This preserves separability 
of the problem of orbital motion when the potential coefficients J 2 t and J 3 are taken into account 
(Reference 3). The inclusion of J 3 is of considerable practical importance, permitting a more 
accurate treatment than that given by perturbation theory. This leads to the computing procedure 
for obtaining the position and velocity coordinates of a drag free satellite from a knowledge of its 
initial conditions. 

STATEMENT OF THE PROBLEM 

If we take r e as the earth’s equatorial radius, and if 

C 2 = r e 2 J 2 (l-TJ3 2 J 2 ' 3 )’ 

and 

s = , 


then the gravitational potential 


V = - fj.(p 2 + c 2 T 7 2 ) 1 (p + T 7 §) 


1 


I 


leads to a separability of the problem of satellite motion. Here, a % + 7 kilometers for the earth 
and the above potential leads to a fit of 


V 


r 



Jn P n (sin O') 


exactly through the third zonal harmonic, about two-thirds of the fourth zonal harmonic, and negligi- 
ble values of order J 2 3 for the higher harmonics. 

If we take p 9 v, and <56 as oblate spheroidal coordinates satisfying the equations 


x + iy - r cos 0 

~ \[p 2 + c 2 )(l 1/2 , 


and 


r sin 6 - prj 


then x, y, and z are the rectangular coordinates of a satellite in a Cartesian frame, with the origin 
at the center of mass of the planet. Also, r , 6, and 4> are the planetocentric distance, latitude, and 
right ascension. 

From Vinti (Reference 4), if a, is the energy of the system, a 3 is the z component of angular 
momentum, and a 2 is the separation constant, then the generalized momenta are given by = a 3 , 


P p - ± (p 2 + c 2 ) 1 F 1/2 ( p ) , 

and 

P,, - ± (l " V 2 )' 1 G I/2 (tj) . 


Here F(p) and G(v) are the quartics : 


F(p) - c 2 a 2 + (p 2 + c 2 ) (- a 2 2 + 2pp + 2a 1 p 2 ) , 


and 


G(77) = " a 2 + (l-7j 2 )(a 2 2 + + 2a t c 2 ?7 2 ) . 
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The Hamilton- Jacobi function W(p, v 9 <t>) is then 


W 



d <fi + 



or 


W - a 3 4> + f ± (p 2 + c 2 ) 1 F I/2 (p) dp + f ± (l - 7? 2 )“ 1 G V2 ( T) ) d ?7 


If fiz and £ 3 ' are constants of the motion, the orbit is then given by 


„ dW 


^2 < 9 a 


and 


_ aw 

^3 " Ja, 


From Vinti's solution of these equations (Reference 3 ), together with the expressions for the gen- 
eralized momenta, we can describe a computational procedure similar to that described in Ref- 
erence 6. 


COMPUTATIONAL PROCEDURE 


Enter the initial conditions x t , y . , z. , x. j y . and z. for a time t . , with the constants p, r e , 

J2 9 J3 < and 

C2 = r e 2 J 2 (l"i J 3 2 J 2 ' 3 ) and 8 = " 5 ' r eJ 2 ' , J 3 > 0 . ( 1 . 1 ) 

Compute 


r . - yx. 2 + y. 2 +z. 2 , r.r. - x.x.+y.y.+z.z., 

i f i J i i ’ ii ii* 7 i>'i ii> 


( 1 . 2 ) 


r? + 2Zj 8 + 6 2 - c 2 j 

n 'S 1 + 


V (r? + 2 z. S + S 2 -c 2 ) 2 


( 1 . 3 ) 




1 



I I II 


II II I II 


2 ( Z i + 5 > 2 

, (the sign of 77. = sign of (z. + b) t 


vt r 


(1.4) 


4 = 2^ (r i + + 


< r i*i * HKj + 22 . 8 + 8 2 - c 2 ) + 2c 2 ( z . + 5)z 
yc r? + 2z. 8 + S 2 - c 2 ) 2 + 4c 2 (z. + S ) 2 


(1.5) 


and 


\ = 


2c 2 77 


, . . . x (r.r. + Sz.)(r 2 + 2z. S + 8 2 _ c 2 ) + 2c 2 (z. + 8 ) z 

[ r i r. + ezj ) + — - - ■ 

y^r 2 + 2 zj S + S 2 - c 2 ) 2 + 4 c 2 (z. + 8) 2 


( 1 . 6 ) 


Compute 


o U 2 


v[Pi + Vi s )(p i + c2r >\ ) 


, 2 \ -1 


(1.7) 


where 


U 7 - x. 2 + y. 2 + z . 2 , 
a 3 = x i ^1 " *i - 


and 


= (1- 


V 


[K 2 


+ c 


2 


) 2 V 


+ a 


2 _ 


[l-7 ? i 2 )(2ai c2 7 7 i 2 


2 m 1 ? ; § )] . 


Then 


and 


_ _ 1 _M 

a o 2 a l e o 



2 a i a 2 2 


1/2 


p 0 = a o( 1 " e o 2 ). 


( 1 . 8 ) 


Prime Constants 

Compute 


X D = -2a. 1 a 2 n~ 2 and X D 4 


( 2 . 1 ) 
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P „ 2 and V = ('sf) • 


K o = TT a "d K 0 J , 
p o 


( 2 . 2 ) 

(2.3) 


[Pi + P 2 ) ' 2P 0 X D - 2 [l-K 0 X D 2 Y D 2 -K 0 2 X D 2 Y D 2 (2 X d 2 -3X d 2 Y d 2 -4 + 8Y d 2 )] . (2.4) 

Pi = Po 2 X d 2 [l + *0 Y d ( x d - 4) - K 0 2 Y„ 2 (12X d j - X D < - 20X D 2 Y„ 2 - 16 + 32Y„ 2 + X D < Y„ 2 )] , (2.5) 




g - 


(P, + P 2 ) 2 ’ 

e = Vl - g , 


, £4> 






-2 = 


2ct j 


2(a 2 2 -a 3 2 ) 


, 1 + 


1 + 


8a, c 2 (a 2 2 -a 3 2 ) 


( c 


2 - Oa r » 2 \ 2 


and 


S = (sin 2 I ) - 


2<Xj c- 


V 


) 2 J 


1/2 


Mutual Constants 

Compute 


p ~ a 


(l - e 2 ) 


- 2ac 2 (ap - c 2 §)(l-S) + 


8a 2 c 2 


S 2 


{l + ip (3S-2)]-S(l-S) 


( ap - c 2 ) (ap-c 2 §) + 4a 2 c 2 S + 5 2 (3ap - 4a 2 - c 2 )§( 1 - §) 

P 


a - 2 A • p o 


(2a)“ 1 (ap''c 2 )A + c 2 , 
(b + ap ~ 2Aa - c 2 ) 


*2' = (^o') 1/2 , 


(jj* (»-«)(»- 17778) 

[l + ( c2 / a o' Po') (1-2S)] 2 


1 

u = 1 + 


U&) 


(1-S) + e , 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 


(3.1) 

(3.2) 

(3.3) 

(3.4) 

(3.5) 
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C 2 

c 2 - 1707 u, 

! c 2 S ^ l 2S „ l c 2 U \ 

c * = PV u (‘-w) • 

(3.6) 

p = 

/ c 2 §uV 1 S - 

l^W) Po' U(1_S) ’ 

(3.7) 

V 0 ' = P + (P 2 

+ S) 1 " , vi = P- (p 2 +§) 1/2 , 

(3.8) 

and 

S' = - Vo vi ■ 

(3.9) 

Using S', we now repeat steps (3.2) through (3.9) to obtain the quantities * 


A, B, a 0 ' , p 0 ' a 2 ' 

, 6 , U, c 2 , Cj , P, 7 7 0 \ Tjj' and S , 

(3.10) 

b | _ “ 2 ^ ’ 

b 2 = B 1/2 , and . 

(3.11) 

Using the formulas of Reference 6 (pages twelve to fifteen) we compute 


Aj = (l - e 2 ) 1/2 

n=2 



and 

A 2 = (l-e»)V*p-l Pn^Rn [( 1 -; 2 ) 1/2 ] . ( 3 . 12 ) 

n=0 / 

1-0 

where P n (b t /b 2 ) is the Legendre Polynomial of degree n, R n (X s ) = X s n P n (x s _1 ) is a ploynomial of 
degree (n/2) in x s 2 , and x s = (l - e 2 ) 1/2 . 

If m is an even integer, compute 




•Equation ( 3 * 9 ) for S* together with the one step iteration (3-10) has now provided for us an accurate value of the element S that was 
originally approximated by T ) Q 2 (Equation 2.10). 
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If m is an odd integer, compute 


^2i+l y ' ( i) 1 " (p 


2i-2n/ b 2 ^ 2n+1 /bj > 

P 2n + 1 ^ > 


Then 


00 

A a = (l _ e2)l /2 p -3 £ Dm R m+2 [{l-e*)i/’] , 


= T U ” e 2 ) 1/2 P~ 3 e (- 2b x b 2 2 p + b 2 4 ) , 


A 12 “ 12 ( 1 " e2 ) 1/2b 2 4e2 P" 3 » 


3 3 


= (l-.»)V.p-ie[b lP -i t ( 3 b*- ba .) p -.-| b|b >( 1+ t)p' 3+ F b 2 4 (4 + 3e2)p-4 

= (l-e*)«/*p-l (abl *- b *)p-«-S e 2 hj b 2 2 p -3 + b* (6e a + e 4 ) P -4 ] , 


A 2 3 " (l-e 2 )" 2 ?- 1 \ (_ b i bj 2 p -3 +V>2 4 p - 4 ) f 


A 2 4 = 256 (l-e 2 )»/ 2 p-5b/e4 , 


A 3 i - (l ~ e2 ) l/2 p' 3 e|^2 +b t p- 1 (3 + f e 2 ) -p" 2 b 2 2 + c 2 ) (4 + 3e 2 )j , 


A 32 - (l" e2 ) 


1/2 !-v - 3 


e 2 3 


T + 4 b iP e "P 






A 34 " " ^ (l-e 2 ) 1/2 p- s e 4 (^b 2 2 +c 2 ) , 


and q = (p 2 +s) l/2 . 




f 

- 


(3.13) 
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Then through third order in j 2 : 

b 2 = 1 - JC.P + (fc 2 + ic 2 )(p 2 + 4q 2 ) + M c 2 2 Q 4 

- I c iC 2 PQ 2 + ^ C 2 C 2 Q 4 + 35§-C 2 3 Q 6 , 


(3.14) 




4 Q S + P 2 - 4 Cj P Q 2 + 4 c 2 P 2 9 2 + (4C 2 + 3Cj 2 ) Q 4 


--3§" c i c 2 PQ 4 + 256 (6C 2 2 + 15C 1 2 C 2 )Q« + ^ C 2 3 Q«(3.15) 


s 3 ° 2 C 2 8 C 1 2 (l6 c i 2c 2 + 8 c 2 2 )( 1 + 2 Q 2 ) 

-i|c 2 3 (i +4 Q 2 + f Q 4 ) +TC,c a P. 

S = r^S’ h, = 4(1 + C 1 -C 2 )- 1/2 , h 2 = 4(1 -Ci-C 2 ) 


- 1/2 


(3.16) 

(3.17) 


and 


- 2 = Q(l-P)- 1 , e 3 = Q( 1 + P) _ 1 , 




2 ttv 1 = ^-2oYJ) 1/2 (a + bj + A t + c 2 A 2 Bj' B 2 ' 1 ) " 1 , 

U' 1/2 A 2 B 2 _1 (a + bj + Aj +c 2 A 2 B/ B^ 1 )' 1 , 


“ a 2 


e' “ a + b where e' < e < 1 


( sy/2 
1 ” B / 


(3.18) 


(3.19) 

(3.20) 

(3.21) 


Here sgn a 3 £o for a direct or retrograde orbit. 

Jacobi Constants 

Compute 


B u = -2PQ + 8 CiQ 3 


(4.D 


8 


(4.2) 


B 12 = - (¥ + ic 2 Q 4 ), 

B 13 = - Cj Q 3 /24 , 


B 14 = C 2 Q 4 / 64 , 

B 21 = -C 2 PQ + ^ C 2 C 2 Q 3 +^-Cj Q 

B 22 = - J2 [(4C 2 + 3C 1 2 )Q 2 +3C 2 2 Q 4 ] 

B 23 ” Te C 1 C 2 Q 3 . 

B 24 = 256 C 2 Q 4 > 


cos <fi. 


yp. 2 + C 2 ^ 1 “ 7 ]? 


sin0 i 


y £ 


l[p} + c 2 y 1 - 17 i 2 


p . 2 + 77 2 C 2 

2 = — 

p. 2 +c 2 


h 2 


p. 2 + TJj* C 2 


cos E. 


sinE. 


M-3)- 

^i h /. (>°i 2 + c2 ) 


ae y/- 2a j' (p. 2 + + B) 


Vi - P 

Sin ^i = Q 

_ ^(p 2 + c 2 T) 2 ) / u 

cos ' Pi Q a 2 y{ i+c lV -c 2 v 2 ) 


(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 

(4.8) 


(4.9) 


(4.10) 


(4.11) 


(4.12) 
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11-111 11 111 11 11 mu mi 1 11 11 11 mu mu iiinii 11 irnin am 



and 


Then 



cos E i “ e 

cos V. = 

1 - e cos E i 

sin v 4 = 

(l - e 2 ) 1/2 sin E a 
1 - e cos E. 

sin nv £ 

for n = 2, 3, 4 

sin n0. 

for n = 2, 4 


cos 3/>. . 


cos e 2 ; 


e 2 + cos 


1 + e 2 cos 


)_ 

(v^r) ’ 

, 5 * 


sin 


(l - e 2 J ) 1/2 sin (■/», + -j) 
('/'i + ?)" 

b-v. P 


1 + e 2 cos 


cos E 3 'j 


e 3 + cos 


$ /r V' 

W-f) 


c 


sinEj'i 


l + e 3 cos(^ 1 - Y) 

S M- ^ 

_ Cot 

(l ~ e 3) U2 Sin _^i_ll) 

l + e 3 cos {&-%) 


(4.13) 


(4.14) 


(4.15) 


E 2i 


E 3i 


Xoi 


Vu 


2 2 Vl + 2S 


E 2i 


E a'i 


2 /1-2C 2 /l>2r 

/3 t = (- 2 a 1 ') _1/2 [bjE; +a(Ej -e sinE ; ) + AjVi + A n sinv. + A 12 sin2v;] 

+ g 2 [b; ■/,. + B n cos i/». + Bjj sin 2p i + B 13 cos 3^ + B u sin 4^]“ 
a L 


(4.16) 


(4.17) 
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^2 a 2 (r 2a iJ 


1/2 




A 2n sinnv. 


+ U 1 / 2 [^B 2 + B 2 j cos ip . 


4 ^ a ' r 

P 3 = 4>. + c 2 a 3 ' (-2a^j' 1/2 A 3 v. + 2^A 3j sinj Vi - U 1/2 |( 1 - S) _ 1/2 [(h 2 + h 2 ) Xq . 

L j=i J 


+ ( h l _h 2 )vii 

] + B 3 ^i- Jc i C 2 Qcos0 i + ^C 2 2 Q 2 sin2A i ) 

.(4.19) 

o 

11 

a 2 B 2 ) ’ 

(4.20) 

and 



l o + e 0 = 277 v i [X 

+ +t>1 +A l) A 2 l ] • 

(4.21) 


Orbit Generator 


Compute 

and 


M s - ( 2"Vjt 

0s = + e„ + 2 71 v 2 t . (5.1) 


By Newton- Raphson iteration we solve M s + E 0 - e' sin (m s +e 0 ) - M s where M s +E 0 - 8; therefore, 


e 


and 8 n = M s initially. 


Compute 


8 . - 


[e n -e‘ s in £ n - M s ] 1 f g n ~ e' sin £„ - Ml 2 |~ e'sin8 n 1 

(l “ e' cos S n ) 2 l-e'cos£ n J j^l - e # cos 8 n J 


(cos£-e) (1-ecos 8) 1 , 


sin v' = (l - e 2 ) 1/2 ( 1 - e cos S)~ 1 sin £ , 5 //t ^ ' O 


(5.3) 


/ - e ' 


11 



; 0 



v o - v ' " M s , 

(5.4) 


(-aa 1 ')-^U-«/*A J B a -‘v 0 . 

(5.5) 

M, = - («+b,)-‘ [ 

( A, + c 2 A 2 B,' B 2 _1 ) v # +^(-2 a,') 1/2 U" 2 B , 2 sin (tf. + ^ )] . 

(5.6) 


M, e * M, 2 sin (M s + E 0 ) 


E 1 = 

cos (M s + E Q ) 2 ^_ e , cos(Ms+Eo) ]3 ’ 

(5.7) 

cos v" 

- [^cos (£ + E t ) “ ej 1^1 - e cos ( £ + Ej )J \ 


sin v" 

= (l-e 2 ) 1/2 [^1 - e cos (8 + Ej )J 1 sin(8 + E 1 ) 

(5.8) 


< 

H 

< 

1 

< 

o 

+ 

3Z 

ll 

< 

1 

< 

(5.9) 

B 22 B 2 -'sin (2/-. + 2^ 0 )+ a 

2 (-2a,')- I/2 U- 1/2 B 2 -‘ [AjV, +Aj,sin(M. + v 0 )+A 2J sin (2M, + 2v 0 )] 



HBj, B' 1 cos (0 S + to) , (5.10) 

M 2 = - (a + b,)-'^, v, +A n sin (M s + v 0 ) + A 12 sin (2M S + 2v 0 ) + (-2a,' ) 1/2 U 1/2 {b,' 0, 


+ B,, cos(0 s +0 O ) 

+ 2B I2 0, cos (20 s + 20 o ) +B ,3 cos (3/> t +30 o ) + B,„ sin (4/>, + 40 o )}j 

- (5.11) 


Mj 

E 2 “ 1 - e' cos (M s + E 0 + Ej ) * 

(5.12) 

cos v'" - 

^cos (£ + E,+E 2 )-eJ|^l-e cos (8 + Ej + E 2 )J 1 , 


sin v'" - 

(l -e 2 ) 1/2 [l ~e cos (8 + E t +E 2 )j _1 sin (£ + E 1 + E 2 ) , 

(5.13) 

v 2 = v "' - ( V 0 + M , +v i) = ( v '"'“ v ")» 

and 

^ = “ B 2 _1 [^21 Sin (^ s + A>) +2B 22 ^cos (2P s +2P 0 ) +B 23 cos (3* b + $P 0 ) + B 24 sin (# s + # 0 )] 

(5.14) 

+ a 2 'u- 1/2 

(-2a, ' )-" 2 B 2 -> [A, v 2 + A 2 , v, cos (M. + v„) + 2A j2 v, cos (2M, + 2v 0 ) 

+ A 23 sin (3M^+3v 0 ) +A 24 sin (4M„ +4v 0 )J . 

(5.15) 
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I 


Then 


E = ? + E t + E 2 , 


V 


M + v „ 


+ V 1 + v 2 


and 


Next 


^ " '/'s + '/’o + ^1 + ^2 


(5.16) 


cos Ej = 


'(** 5 ) 


(L 6 b C o >, ft t 


1 + e 2 cos 


sin E 2 


(l - e 2) 1/2 sin ^+j) 

1 + e 2 cos (xf> + y j 


cos E/ 


e 3 + cos 




1 + e 3 cos 


(*-?) ' 


sinE 3 


l 1 " e 3 2 )‘ /2 Sin (’/’* T) 


1 + e 3 cos 


(♦-?) 


(5.17) 


Vo 


V E 3 

+ 


2 yi - 2c 2 yi + 2 ? * 


Vi 


E 2 ' 


2 y 1 - 2 s 2 yr+2T 


(5.18) 


P _ a( 1 - e cos E) , 17 - P + Q sin < p 
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4> - P 3 ~ C * a 3 (~ 20L 1*) ~ 1/2 A 3 v+ /* A 3j Sin jv + ^7 ° 1/2 {(l~ S )~ 1/2 [( h l + h 2 ) Y 0 


+ ( h l ~ h 2 ) Yi] +B 3 C 1 C 2 Qcos0 + 33 c 2 2 Q 2 sin 20 j- , (5.19) 


, 2 _ P^_ C 2 

^ o 2 +n 2 


p 2 + 7] 2 C 2 


V - (p 2 + c 2 )( l-v 2 ) , 


(5.20) 


ae y- 2 aj' (p 2 + Ap + B) 

A = h V(,..c^ *‘" E • 


Qa 2 # |/l + C x 77 - C 2 T) 2 cos \p 

t/u ( p 2 + C 2 7J 2 ) 


(5.21) 


X = /(p 2 + C 2 ) ( 1 - tj 2 ) COS 0 , 


Y = y(p 2 + c 2 ) ( 1 -T 7 2 ) sin 0 , 


PV ~ 5 . 


( PP 3?Z7_\ _ _ 

\p 2 + c 2 1 -rj 2 j h 

/ p'p + 

\p 2 + c 2 1 — T 7 2 / h 


*-«3* 


X-a 3 ' 


^ ’ 


Z = p-q + Tjp . 


(5.22) 
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CONCLUSION 


The accuracy of the orbit itself as a solution for the given potential is carried out through 
terms of the third order in J 2 , the coefficient of the second zonal harmonic. Its accuracy however, 
and thus that of the secular terms, may be increased at will. Periodic terms are carried through 
the second order, but their accuracy may also be increased. In order to carry this kind of accuracy 
most perturbation methods presently in use would become far too cumbersome and impractical for 
orbit computational purposes. 

An advantage of accounting for J 3 in this way is the absence of small denominators in e or 
sin I that occur in certain perturbation theories. Thus, one can easily compute polar orbits and 
circular equatorial orbits. 

This program, similar in structure to the previous Vinti accurate intermediary orbit requires 
a relatively small number of storage locations throughout the entire computing procedure. Recent 
tests on the IBM 7094 have indicated that this latter Vinti method enjoys a rather rapid computa- 
tional speed and it is indicated through preliminary comparisons that it is faster than other methods 
which use perturbation techniques to compute orbits. Electronic timing on the Theoretical Division 
computer shows that this program can compute 3000 minute vector points (time, x, y, z, x, y, and z, 
3000 times per minute of IBM 7094 computer operation) with simultaneous production of BCD tape 
output (Reference 6). 

The residual fourth harmonic term (Reference 5) has been programmed and is presently being 
tested as is the J 2 t tesseral harmonic term. In addition, the orbit differential correction program 
written for the Vinti accurate intermediary orbit theory (Reference 7), will need only slight modifi- 
cation to account for the element S which replaces = sin I in Vinti T s new theory. Consequently 
it is expected to go just as rapidly, producing a set of mean elements of even greater accuracy with 
which to predict satellite orbits over a longer interval of time. In addition, comparison with a very 
accurate, double precision numerical integration for polar and near polar (89° inclination) and 
equatorial cases over a five day period indicate there is no loss of accuracy from the general or 
intermediate inclination cases. The agreement between the Vinti program and the numerical inte- 
gration was always six to eight decimal places. The Vinti method computed in single precision or 
eight decimal places only. 

It has been shown (Reference 7) that the orbit differential correction program for the Vinti 
accurate intermediary orbit theory converges rapidly and with great accuracy. This latter program, 
tested on the Relay II Satellite with an eccentricity of 0.23597617, and Satellite ANNA with an ec- 
centricity of 0.00671710 gave the following results: 

Using radio direction cosine observation data, and for seventy equations of condition of the 
Relay II Satellite extending over a five hour arc following injection into orbit, the program con- 
verged on the third iteration to a standard deviation of fit criterion of 2.7 x 10 ~ 3 within thirty 
seconds. Using Smithsonian Astrophysical Observatory optical data from the Anna Satellite ex- 
tending over an arc of seventy-five hours following injection, and with twenty equations of condition, 
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the program converged on the second iteration to a standard deviation of fit of 0.2 x 1CT 3 . For 
seven equations of condition of Satellite ANNA covering the first forty-five hours following injec- 
tion, the Vinti orbit differential correction program converged on the second iteration to a standard 
deviation of fit criterion of 0.04 x 10 ~ 3 . All tests were conducted using the IBM 7094 Mod. I elec- 
tronic digital computer. The orbit differential correction program utilizing the new Vinti equations 
of motion, now carries the additional accuracy inherent in the new method, and should merit even 
tighter tolerances for the standard deviation of fit. 
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